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Enrico Carlon 1 and Marco Baiesi 2 

'interdisciplinary Research Institute c/o IEMN, Cite Scientifique BP 69, F-59652 Villeneuve d'Ascq, France 
2 INFM Dipartimento di Fisica, via Marzolo 8, 1-35131 Padova, Italy 
(Dated: February 2, 2008) 

We analyze the unbinding transition for a two dimensional lattice polymer in which the constituent 
strands are mutually avoiding random walks. At low temperatures the strands are bound and form 
a single self-avoiding walk. We show that unbinding in this model is a strong first order transition. 
The entropic exponents associated to denaturated loops and end-segments distributions show sharp 
differences at the transition point and in the high temperature phase. Their values can be deduced 
from some exact arguments relying on a conformal mapping of copolymer networks into a fluctuating 
geometry, i.e. in the presence of quantum gravity. An excellent agreement between analytical and 
numerical estimates is observed for all cases analized. 
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I. INTRODUCTION 

The unbinding transition from a low temperature 
double stranded polymer to a high temperature sin- 
gle stranded phase has been the subject of recent at- 
tention in the context of studies of DNA denaturation 
|H Q EJ El S Two main approaches have been used to 
model this unbinding. The first one relies on the use of 
directed polymers, where only the transversal coordinate 
measuring the distance between homologous base pairs is 
considered In a second approach, one considers the 

polymer as being composed of an alternating sequence of 
double stranded segments and denaturated loops 0, @ • 
The statistical weights assigned to loops and segments 
can be estimated using concepts of homopolymers and 
self-avoiding walks (SAW) statistics. 

Traditionally, in the latter class of models, the statis- 
tical weight of a loop was approximated as the number 
of configurations for a closed SAW, neglecting any ex- 
cluded volume interaction with the rest of the chain . 
More recently, statistical mechanical ideas based on the 
theory of polymer networks [lflj , were used to take into 
account the excluded volume effects in an approximated 
way 0- Although quite simple, this analysis, which was 
carried out analytically, captures quantitatively very well 
the asymptotic form of the loop partition function, as a 
series of numerical investigations on two and three di- 
mensional lattice models have shown 

HGH]. Des P ite the 
general good agreement between analytical predictions 

and simulations, quite small but systematic deviations 
between the two were found . 

The aim of this paper is to investigate further on these 
issues for other types of models of polymer unbinding. 
We focus on a two dimensional lattice model for which 
we obtain a series of analytical predictions based on ex- 
act results from conformal invariance |12j . In this model, 
the two constituent strands are two random walks (RWs) 
with an attractive interaction. While we relax the ex- 
cluded volume interactions within each strand, mutual 
avoidance, i.e. the non-overlapping condition between 
the strands is preserved. This makes the bound dou- 



ble stranded state to behave as a SAW, assigning to it 
a quite different physics compared to that of loops. We 
show that in this model unbinding is a very strong first 
order transition. 

As is well-known from polymer physics, the partition 
function of a closed SAW of total length I assumes the 
following asymptotic form ^| 



Z(l) ~ (J, I 



l]— Co 



(1) 



where /i is a geometric factor and Co a universal expo- 
nent which equals cq « 1.76 [j| in three dimensions. It 
has been shown 0] that a loop attached to two long seg- 
ments or loops (see Fig.^) has still a partition function of 
the form of Eq. , but with different exponents. For in- 
stance in three dimensions one finds for a loop embedded 
between two long segments Q c ss w 2.1 and c\\ w 2.2 for 
a loop embedded between two long loops (here we used 
the subscripts s or I to indicate neighboring segments or 
loops). An increase of c for embedded loops is caused 
by the tendency of the loop to become more "localized" 
due to excluded volume interactions with the rest of the 
chain. For c > 2 the transition becomes first order 0, [j| . 

In this paper we calculate analytically the values of 
the exponents cn and c ss , as well as other type of en- 
tropic exponents, for the model of mutually avoiding ran- 
dom walks. The calculation relies on conformal map ping 
of copolymer networks into a fluctuating geometry [12j , 




FIG. 1: Hierarchy of exponents describing the loop entropy 
for an isolated loop (co), a loop embedded between two long 
segments (c as ) and a loop embedded between two long loops 
(co). 
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FIG. 2: Snapshot of a configuration for the model studied 
in this paper generated by the Pruned Enriched Rosenbluth 
Method (PERM). The two walks are indicated as solid and 
dashed lines. When they are bound only the solid line is 
visible. Inset: Blow up of a small area of double stranded 
phase with two embedded loops. The loops are oriented, i.e. 
following the direction of one of the strands, the other strands 
is always found on the same side. 

i.e. in the presence of quantum gravity. This theory, 
which was recently applied to the calculation of multi- 
fractal spectra of harmonic measures 0J > enables one to 
obtain exact entropic exponents of networks composed 
of arbitrary mixtures of random and self-avoiding walks 
[l2| . We also performed a series of Monte Carlo simula- 
tions in order to determine the numerical values for the 
exponents. In all cases analyzed the numerical and an- 
alytical results are in excellent agreement. This is due, 
as we will discuss below, to the strong first order nature 
of the unbinding transition, which makes the model an 
ideal testground where copolymer network theories can 
be applied. 

Part of the results presented here, have been discussed 
in concise form in Ref. 14] . In this paper we present 
the results of large scale numerical calculations, which 
we extend to other quantities not considered previously, 
and present a full account of the analytical results. The 
present model has also been investigated by means a con- 
tinuum approach £| , in which mutual avoidance between 
the two strands has been approximated by an effective 
long-range interaction, an approach which predicts a first 
order unbinding transition. 



II. THE MODEL 

We consider two random walks of length N on a 
square lattice described by the vectors f\{k) and r<2(fc) 
with k = 0, 1 ... N. The walks have common origin 
ri(0) = ^2(0) and are not allowed to overlap except at 
homologous sites, i.e. ri(i) = r*2(j) is possible only if 
i = j. Whenever such a contact is realized the sys- 




FIG. 3: Example of a star copolymer formed by two random 
walks and a self-avoiding walk, where all three of them are mu- 
tually avoiding and have average size equal to R. The exact 
entropic exponent for any star copolymer with an arbitrary 
number of random and self-avoiding walks can be calculated 
thanks to a mapping onto a fluctuating geometry (see Sec. 

urn . 



tern gains an energy e = — 1. At very low temperatures 
the two walks are fully bound and form a self- avoiding 
walk, since, as the walks are mutually avoiding, a bound 
site fi(i) = ^2(2) cannot overlap any other sites on both 
strands. At higher temperatures unbinding starts from 
the unconstrained edge i = j = N and loops proliferate 
along the chain. Figure [21 shows a snapshot of a con- 
figuration generated by the Pruned Enriched Rosenbluth 
Method (PERM) where the two strands are rep- 

resented by solid and dashed lines, respectively. Notice 
that the double-stranded part has a characteristic self- 
avoiding walk behavior and the unbound single strands 
behave as mutually avoiding random walks. The inset of 
Fig. [2] shows a blow-up of part of the double stranded 
chain with two short loops of few lattice spacings of 
length. The model is constructed such that following 
the two strands ri(fc) and r2(k) from the origin k = to 
k = N one finds one of the two strands (say ri(k)) when 
unbound always at the left side of the other strand, i.e. 
the loops are oriented. This choice will allow us to re- 
strict the type of diagrams considered in the continuum 
polymer network description of the model. 



III. CONFORMAL MAPPING ONTO TWO 
DIMENSIONAL QUANTUM GRAVITY 

A. Star copolymers 

Duplantier formulated an elegant theory which al- 
lows one to compute exact entropic exponents for star 
copolymers, an example of which is shown in Fig. [3J A 
star copolymer is formed by an arbitrary mixture of ran- 
dom and self-avoiding walks joined at a common origin, 
which can be all mutually avoiding (as in the example of 
Fig. |3| or partially transparent to each other. 
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One is typically interested in the grand-canonical par- 
tition function Z R (S), for a star S formed by /i self- 
avoiding walks (SAWs) and f% random walks (RWs), all 
with average size R from the origin of the star, R denot- 
ing the end-to-end distance for each walk. Having the 
same size, SAWs and RWs are characterized by different 
lengths, scaling as ~ R 1 ^, with v = 3/4 for a SAW and 
v = 1/2 for a RW. Fugacities per unit of step length are 
associated to each type of walk. By tuning their values 
appropriately both walks become critical. In this limit 
the partition function scales asymptotically for large R 
as 12] 



Z(S) ~ 



(2) 



where r](S) is the scaling exponent associated to the sin- 
gularity at the center of the star and 7/2 the entropic expo- 
nent associated to an isolated SAW (772 = —11/24 in two 
dimensions, while the corresponding exponent for RWs 
is zero) . We have followed here the notation of Ref . ^(| , 
which is slightly different from that of the original work 
of Duplantier [lj| , but more suitable for a generalization 
to networks. The star exponent can be written in the 
following form: 



r)(S) = -2A(5) + /1 



2 ' 



(3) 



with A(S) the conformal scaling dimension and ry^ = 
5/24 the correlation length exponent |lq . 

Here we review briefly the main formulas leading to 
the exact value of the exponent A(S) for an arbitray 
star copolymer S. Details of the derivation can be found 
in the Refs. [T2I Il7| . The main idea is to map the star 
copolymer from the planar euclidean geometry onto a 
two dimensional random lattice, i.e. in the presence of 
quantum gravity. 

Besides the bulk scaling dimension A(S) of a star S 
we will also consider star copolymers confined in a half- 
plane with the origin near the boundary of the plane, 
which defines the surface scaling dimension A(S). We 
use the notations A qg (S") and A qs (S') for bulk and sur- 
face conformal dimensions in the random lattice (here qg 
stands for quantum gravity). Given two walks A and B, 
with a common origin we indicate (as in Ref. [l^) with 
the symbol A V B a star configuration where the two 
walks are allowed to overlap each other and with A A B 
a configuration where A and B are non-overlapping. 

The main result of the theory is that in the fluctuat- 
ing geometry the surface conformal dimensions for two 
mutually avoiding walks are additive |l2l | i.e.: 



A qs (A A B) = A^(A) + A qg (5). 



(4) 



For transparent walks in the plane, due to the trivial 
factorization of their partition functions one has: 



A(A V B) = A{A) + A{B). 



(5) 



We point out to a sort of duality between Eqs. (JIJ and 
(J5J. Eq. states that non-intersecting walks, become 




FIG. 4: Star copolymers considered in this paper. Thin lines 
denote random walks, while thick lines denote self-avoiding 
walks. Solid and dashed thin lines are allowed to overlap 
with thin lines of the same species, but not allowed to over- 
lap thin lines of the other species (eg. overlap between solid 
and dashed thin lines is not allowed). The conformal scaling 
dimensions for the examples in the figure are A(S a ) = 5/8, 
A(S b ) = 39/32 and A(S C ) = 35/24 (see text). 



transparent to each other when placed in a fluctuating 
geometry. The Eqs. (@J and JSJ) can be generalized to 
any nested star copolymer structure (see Ref. [l2|). 

Conformal invariance relates the bulk and surface scal- 
ing dimensions in the two geometries as |T^: 



A = [/(A qg ), 
A = C/(A qg ), 



(6) 
(7) 



where U(x) 
i(V24¥^T 



f (1 + 2x) and the inverse U 1 (x) = 
1). The previous equations follow from 



the work of Knizhnik, Polyakov and Zamolodchikov , 
who established the existence of a relation between the 
conformal dimensions of scaling operators in the plane 
and those in the presence of gravity. Finally the relation 



A qs = 2A° 



1 

2' 



(8) 



connects surface and bulk conformal dimensions in the 
fluctuating geometry and can be derived from some fac- 
torization properties of star partition functions under 
quantum gravity |19|. 

For random and self-avoiding walks the conformal di- 
mensions arc 121: 



*RW 



1. 



A qg 



1, A 



1 



RW 



A qs 

^RW 



1 



(9) 



Asaw = I Af m = \, A SAW = |, Af AW = i (10) 



As a practical example of the use of the above formulas 
we calculate the conformal dimensions of the three star 
copolymers shown in Fig.0J which are the configurations 
relevant for the model discussed in this paper. 

Let us consider first the star copolymer composed by 
two mutually avoiding random walks (i.e. S a of Fig. [3J. 
The equation J2J states that the surface conformal di- 
mension in the fluctuating geometry for two mutually 
avoiding walks is additive, thus for S a it is twice as large 
as that of a random walk 



A qg (S a ) = 2 A 



RW 



(11) 
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From Eq. © one finds therefore for the bulk dimension 
in the fluctuating geometry: A qg (S a ) = 3/4. The final 
step to obtain the bulk dimension in the plane is to use 
the conformal mapping [Eq. ©] which yields 

A(S a ) = C/(3/4) = |. (12) 

The generalization to k mutually avoiding random walks 
is given in the Appendix ^] The above derivation of 
A(S a ) illustrates the general strategy of the calculation: 
the conformal dimension is first calculated in the fluctu- 
ating geometry, where mutual avoidance is easy to im- 
plement [see Eq. @] and then obtained for the planar 
geometry from Eqs. © and J7J). 

For the star formed by two random walks and a self- 
avoiding walk all of them avoiding each other (Sb of 
Fig. EI, one proceeds along the same lines as done for S a . 
First, the additivity of surface scaling dimensions in the 
fluctuating geometry [Eq. ©] implies that A qg (S'h) = 
2A^ + A^ w = 11/4. Therefore, for the bulk dimen- 
sion one finds [Eq. ©]: A^(S b ) = 9/8. Finally, the 
conformal mapping [Eq. ©] yields: 

3Q 

A(S b ) = 17(9/8) = -. (13) 

For the star copolymer S c of Fig. 0] the calculation is 
slightly different. Both solid and dashed walks are trans- 
parent to each other, therefore we first need to calculate 
the scaling dimension for the sub-star composed either 
by solid or by dashed lines only. Two transparent ran- 
dom walks (RWi and RW2) in the planar geometry have 
additive surface conformal dimension [Eq. ©] thus 

A(RWi V RW 2 ) = 2A RW = 2. (14) 

By inversion of Eq. J7J, we obtain the corresponding 
surface conformal dimension in the fluctuating geometry: 
A<K(RWiVRW 2 ) = U~ 1 (2) = 3/2. Now, as in the fluctu- 
ating geometry the two scaling dimensions of the mutu- 
ally avoiding dashed and solid sub-stars are additive one 
finds from Eq. ©: A^(S C ) = 2A q «(RWi V RW 2 ) = 3. 
Eq. © yields for the bulk: A q s(S c ) = 5/4. The fi- 
nal step is the mapping back into the planar geometry 
[Eq.©] which yields: 

35 

A(S C ) = C/(5/4) = -. (15) 



B. Networks with arbitrary topology 

We consider now exponents associated to networks 
of arbitrary topology, formed by mixtures of RWs and 
SAWs connected to each other in such a way to form 
loops and dangling ends (examples are shown in Fig. [S] 
and©. The partition function for a network G contain- 
ing fi SAWs, in the limit when all walks are critical takes 




FIG. 5: Examples of loop configurations relevant for the 
model, (a) Isolated loop formed by two mutually avoid- 
ing random walks, (b) loop embedded between long SAWs 
(R ^> r), (c) Loop embedded between two other long loops. 
At the edges the first loop formed at the origin of the two 
walks can be attached to a long SAW (d) or to another long 
loop (e). 

the form 0] 

Z(G) ~ BPa-hm^ ( 16 ) 

where the universal exponent t\q depends on the topology 
of G as follows [13 

VG = -d£ + Y,nsv(S) (17) 

s 

where d is the dimensionality of the system (d = 2 in 
this paper), C the number of independent loops and the 
sum is extended to all constituent vertices forming the 
network, each contributing for a factor 77(5) (the star ex- 
ponent defined by Eqs. J3J and ©, and ng is the degen- 
eracy of the vertex S. Using the results of the preceding 
section and recalling that for a vertex S with f\ outgoing 
SAWs one has r)(S) = -2A(S) + [Eq. ©] we can 

now calculate the network exponents for the examples in 
Fig. [S] and which arc those relevant for the unbinding 
transition considered in this paper. 

For the isolated loop of Fig. ©a) one has C = 1 and 
there are two vertices S a , as defined in Fig. 21 As there 
are no SAWs j\ = in Eqs. ljlU|l and © therefore Eqs. 
ijnj and (O imply that 

7 ?G = -2 - 27 7 (S Q ) = -2-4A(S' a ) = -2-^ (18) 
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Usually we are interested in the scaling as function of the 
loop length, and not of its radius of gyration, therefore 
the partition function is: 



Z 



loop 



v 



(19) 



with v — 1/2 for random walks and where we have in- 
troduced the entropic exponent for an isolated loop, in 
analogy as what discussed in the Introduction. We have 
co = 2 + 1/4 = 2.25. Differently from the case of self- 
avoiding loops (see Introduction) in the present model 
already at the level of an isolated loop Co > 2, implying 
that the first order character of the transition is rather 
strong, as remarked in Ref. [Til ]. 

Equation l|16[> can be generalized to the case where the 
network is formed by walks of different sizes, say R and 
r. In this case the network partition function becomes: 



Z(G) ~ jpo-hmf (Z.^ 



(20) 



with / a scaling function. In all the other cases (Fig.EJb- 
e)) we will consider the limit in which r, the size of the 
loop, is much smaller than R the size of the walks or 
loops attached to it. 

We consider now the network of Fig. |SJb) which con- 
tains one loop C = 1 and two vertices Sb (as in Fig. |3J| 
therefore: 



Vg - 2?72 = -2 



-2A(S b ) + ?±- m 



(21) 



In the limit r -C R one should recover the partition func- 
tion of a single SAW (~ R~ n2 ), which implies that the 
scaling function of Eq. 1)20(1 behaves as 



fix) 



„-2-4A(S' i ,)+ J ;*-r )2 



for 







(22) 



In this limit the partition function of the network be- 
comes: 



Z(G) ~ ^-»)2 r _2_4A ^ s ' 6 ' +?, * _I ' 2 



(23) 




FIG. 6: Configurations of end-segments linked to a long SAW 
(a) and to a loop formed by two mutually avoiding random 
walks (b). 



and thus factorizes as Z(G) ~ -ZsAW-^ioop, i.e. in a con- 
tribution from a long SAW and of a loop. The latter 
expressed in terms of its total length I reads: 



7 loop 



|-v[2+4A(S 6 )-^+»7a] ^ j- 



(24) 



Using the numerical values of the exponents give n in the 
preceding section we find: c BS = 3 + 5/32 |21|. Note 
that the derivation of Eq. (|24[1 is similar to that for the 
unbinding of SAWs reported in Ref. Q ■ 

The next example is the configuration of Fig. |3Jc), i.e. 
a loop confined between two other long loops. As the 
calculation of its partition function follows closely that 
of a loop confined between two SAWs we report here the 
final result: 



7 loop 



-i/[2+2A(S c )] 



(25) 



and thus from Eq. l(T5|) we find: en = 2 + 11/24. Notice 
that cn< c ss , contrary to what happens for self-avoiding 
walks 2] . This follows from the random walks character 
for isolated strands: As solid and dashed lines in Fig.(5fc) 
are allowed to overlap themselves a loop bounded by two 
loops is "less localized" (smaller c) than a loop bounded 
by two SAWs. 

We will also be interested in the statistical properties 
of the first loop formed at the common origin of the two 
walks, therefore we consider also the case of a loop bound 
to the rest of the polymer only on one edge (see Fig. EI d- 
e)). Also in this case we only report the final results 
as the calculation follows the example above. For the 
partition function of a loop bound to a long segment 
(Fig. Eld)) we find 



v loop 



^-i/[2+2A(S a )+2A(Si,)-^/2] ^ j- 



(26) 



with c s = 2 + 19/24 while for a loop attached to a long 
loop (Fig. Ete)) we find: 



^loop 



£-i/[2+2A(S c )] ^ i~ C1 



(27) 



with c\ = 2 + 11/24. Notice that, curiously, c\ = c\\ a rela- 
tion which is not valid only for this particular model, but 
it is quite general for all polymer unbinding transitions, 
wether the constituents strands are SAWs or mutually 
avoiding RWs. 

Finally besides loops it is also interesting to consider 
end-segments distributions, i.e the length of the single 
strands at the free end of the polymer (see Fig. EJ, as 
done for the SAWs in Ref. (2^ • Again, as the calculation 
is very similar to those reported in this section we only 
give the final results. For end segments each of length n 
bounded to a SAW of size R, in the limit R^$> r, with r 
the size of the end-segments (r ~ n 1 / 2 ) we find: 



^cnd 



-i/[2A(S b ) -7)^/2] 



(28) 



where we have introduced a new exponent equal to 7 S 
7/6. 
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TABLE I: Summary of the exact loop and end-segments en- 
tropic exponents for the configurations of Fig. |S]and|n] 



Loops 




End sej 


;ments 


CO 


2 + 1/4 (»2.25) 


7o 


5/8 


0.62) 


('ss 


3 + 5/32 («3.16) 


7s 


7/6 (ft 


1.17) 


Cll 


2 + 11/24 (« 2.46) 


7i 


5/6 (ft 


0.83) 


c B 


2 + 19/24 (« 2.79) 








Cl 


2 + 11/24 (« 2.46) 









Similarly for the configuration of Fig. Etg) we find 

Z cnd ~ n -[2A(S c )-2A(S a )] „ n - 7! (2Q) 

with 7; =5/6. Notice that analogously as what we have 
found for the loops, also the end-segments bound to a 
long loop are less localized than those bound to a SAW 
(7/ < 7s)- 

The results obtained in this Section are summarized 
in Table [I] to which we have also added 70 the exponent 
associated to isolated end-segments. This exponent is 
associated to the configuration of Fig. Ufa), therefore: 
70 = 2uA(S a ) = 5/8. 

IV. NUMERICAL RESULTS 

Polymer configurations in which each strand is of 
length N = 1280 were generated by the Pruned En- 
riched Rosenbluth Method (PERM) which is described 
in Ref. ^^|. We performed first three runs at fixed tem- 



peratures around the critical point and used the multi- 
hystogram method [23| to interpolate results to arbitrary 
temperatures in the transition region. A precise estimate 
of the transition point was obtained by the analysis of 
the specific heat maximum per unit of length C max (A^), 
which is expected to scale as function of the chains length 
N as [3 

C max (N) ~ N 2 ^ 1 (30) 

a relation which defines the crossover exponent cj>. Figure 
0a) shows a plot of the peak heights C maK (N) as a func- 
tion of the chain length N on a log-log scale. A linear 
fit of the data yields (j> — 1.000(2) in excellent agreement 
with a first order transition, for which one expects = 1 
fl3| . The sharp determination of is a signature of a 
rather strong first order character of the transition. It 
is interesting to point out that in the case of unbinding 
of self-avoiding walks, although the transitions is known 
to be of first order type, it is difficult to extrapolate an 
exponent <fi which is consistent with <fi = 1 [3J, |5j (at least 
in three dimensions). 

The peak position is expected to scale as 

/WA0 ~ Pc + ^ (31) 

with A a constant. The first order character (0=1) 
is confirmed by the scaling of the peak positions, as il- 
lustrated in Fig. 0b) which shows a plot of Pm ax {N) 
vs. 1/N ((3 is the inverse temperature). We performed 
a series of iterated linear fits of f3 max (N) vs. 1/N and 
obtained the following estimate of the transition point 
P c = 1.7617(1). 
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FIG. 7: (a) Plot of the specific heat peaks heights C max as 
a function of the chain length N on a log-log scale. From a 
linear fit and by Eq. (1301 we find as estimate of the crossover 
exponent <f> — 1.000(2). (b) Plot of the inverse temperature 
position of the peaks as function of 1/N. 



A. Behavior at the transition point (5 = j3 c 

We focus first on the behavior at the transition point 
j3 = [3 C = 1.7617. The probability distribution of finding 
a loop of length / is expected to decay as a function of I 

as 5] 



p(i) ~ r 



(32) 



from which the exponent c can be calculated. Figurc^a) 
shows a log-log plot of P(l) versus I for N = 1280 at j3 c . 
As a comparison we plot, as straight lines, the slopes cor- 
responding to the analytical estimates of the exponents 
c ss and en, f° r a 1°°P embedded between two segments 
and two loops, respectively, calculated in the preceding 
section. Notice the excellent agreement of the numerical 
results with the decay exponent c ss , which indicates that 
considering each loop as simply bounded by pure SAWs 
approximates extremely well the polymer configuration 
at the transition point. As P(l) decays rather fast in I 
(~ P c with c w 3.2), rather long runs are needed in order 
to obtain an accurate statistics. 

We have also analyzed the probability distribution of 
the first loop Pq(1), which is formed at the common origin 




FIG. 8: Probability distribution of loop lengths at the esti- 
mated transition point f) = (3 C — 1.7617 (a) and in the high 
temperature phase /3 = 1.6 (b). As references we plot the 
exact exponents for a loop embedded between two segments 
c ss = 3 + 5/32 and a loop embedded between two long loops 
c\\ — 2 + 11/24, calculated in the preceding section. 



FIG. 9: Probability distribution of the first loop as a function 
of its length in a log- log scale at the estimated transition point 
(3 = f3 c — 1.7617 (a) and in the high temperature phase j3 = 
1.6 (b). The slopes corresponding to the analytical estimates 
for the exponents c s and ci for a loop bounded by a SAW 
segment and by another long loops are also shown. 



of the two strands. Here we have considered only loops 
originated at the very edge for which the first monomer is 
unbound ^ ^(l)). A plot of log\Po(0 as a function 

of log I at (3 C is shown in Fig.[5Ja). The statistics is poorer 
compared to the total loop distribution in Fig. OJa), as 
in most of the configurations the first monomer is bound, 
therefore loops at the common edge of the two strands 
are quite rare. Despite that, the agreement with the 
expected exponent c s = 2 + 19/24 is very good. 

Finally we considered the end-segments distribution 
which is shown in Fig. HOf a). Here P e {n) is the prob- 
ability of having a configuration in which the last n — 1 
monomers are unbound while f\{n) = ^(rt) . Once 
again at (3 C we note a good agreement with a decay 
P c {n) r*> n~ ls with an exponent j s = 7/6. 



B. Behavior in the high temperature phase < f3 c 

We have repeated the same type of calculation of loops 
and end-segments statistics also in the high temperature 
region (3 < (3 C . We expect a power-law distribution of 
loop lengths also at high temperatures 

Figure IHIb) shows a plot of the loop probability dis- 
tribution at f3 = 1.6 for N = 1280. Indeed there is a 
clear power-law decay governed by an exponent which is 
in excellent agreement with c\\ — 2.46 as calculated in the 
preceding section. Typically at high temperatures a loop 
is more likely to be bound by neighboring loops, rather 
than double stranded segments, as contacts between the 
two strands are rare, which explains the observed expo- 
nent. 

Analogously, also for the first loop length distribution 




FIG. 10: Probability distribution for the end segments at the 
estimated transition point (3 — /3 C = 1.7617 (a) and in the high 
temperature region f3 = 1.6 (b). The analytical estimates 7 S 
and 7; for the slopes for end segments attached to a SAW and 
to a loop respectively, are also shown. 



Po(l), we find a decay exponent in very good agreement 
with ci (see Fig.^b)). Notice that, differently from the 
[3 — [3 C case, here the distribution is rather smooth, as, 
in the high temperature region, it is more likely to find 
configuration where a loop forms at the origin. 

An analogous very good agreement with the exponent 
7/ has also been found for the decay of the end segments 
distribution P e (n), as illustrated in Fig. Iior bl. 
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V. DISCUSSION 

In this paper we have studied the unbinding transi- 
tion for a two-dimensional lattice polymer composed of 
two strands which are random walks mutually avoiding 
each other. Considering all the polymer unbinding mod- 
els studied so far in the literature in two and three dimen- 
sions, the present model is that with the strongest first 
order transition (higher loop exponent c ~ 3.2). The 
sharp first order behavior can be inferred from the deter- 
mination of the crossover exponent <j>, which is in excel- 
lent agreement with the first order value <f> — 1, already 
clearly observed for rather short chains (N w 100). This 
can be compared with, for instance, the unbinding for 
SAWs in three dimensions for which c ss 2.1 (weak first 
order, just above the threshold c = 2) and numerical 
4> ~ 0.9, in itself not fully consistent with a first order 
transition Q. 

We showed that numerical estimates of entropic expo- 
nents for loops and end segments are in excellent agree- 
ment with analytical results, obtained from a recent the- 
ory based on mapping of star copolymers into a fluctu- 
ating geometry [12^. Such good agreement is not sur- 
prising in the high temperature phase where the loops 
are most likely bound by other loops as contacts be- 
tween the strands are very rare. At the transition point 
one expects that typical configurations are composed by 
a double stranded polymer "dressed" with loops of all 
sizes. Notice however that as the loop exponent is rather 
large (c w 3.2), the statistical weight of long loops is sup- 
pressed. In particular c > 3 implies that the two first 
moments (I) and (I 2 ) are finite. In this case, as loops 
are typically very small (see also Fig. |2J, therefore ne- 
glecting totally their effect is still a very good approxi- 
mation. Thus the strong first order unbinding represents 
an ideal case where polymer networks calculations work 
extremely well. 

As a counterexample we mention the case of a two 
dimensional unbinding transition studied recently in di- 
block copolymers [24( ■ A diblock copolymer is composed 
by two homogeneous branches of A and B monomers 
joined at a common origin. In the model of Ref. [24| A 
and B are both self- and mutually avoiding. An attractive 
interaction between A and B induces a " zipping" transi- 




FIG. 11: Zipping transition for a two dimensional diblock 
copolymer. Both strands A and B are mutually and self- 
avoiding. Loops in the zipped phase may have different 
lengths in the two strands. 



tion by lowering the temperature where the two strands 
are bound (see Fig. 111(1 . Differently from the model stud- 
ied here and in models of DNA denaturation 0, 0, El 13 
any "monomer" in A can bind to any "monomer" in B, 
therefore loops of total length I may have different lengths 
along the two strands, i.e. I = I a + Ib with Ia ^ Ib- Due 
to this freedom the loop entropic exponents for the di- 
block copolymer are given by: 

c ss = 4 SAW > - 1 « 1.42 (33) 
cu = 4 SAW) - 1 « 1.64 (34) 

where we have used the two dimensional SAW exponents 
0: c f s m) = 2.42 and Cl ( f AW) = 2.64. Numerical re- 
sults Q shows that the zipping transition is continuous 
and the specific heat exponent is in excellent agreement 
with (j) = 9/16, conjectured to be an exact value. From 
scaling arguments for a continuous transition Q one has 
c = 1 + cj) — 25/16 ~ 1.56, which is in between the values 
of Eqs. (|33|) and (|34p. clearly distinct from both. At the 
zipping transition therefore the polymer network theory 
does not reproduce the numerical value for the loops en- 
tropic exponent as accurately as for the model studied 
here. Notice that however, also in this case the numeri- 
cally determined c satisfies the relation c ss < c < en, as 
expected. 
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APPENDIX A: STAR COPOLYMER MADE OF k 
MUTUALLY AVOIDING RANDOM WALKS 

We generalize here the calculation leading to Eq. I|12|) 
to the case of consider a star copolymer Sk made of k 
mutually avoiding random walks. In this case the surface 
conformal dimension in the fluctuating geometry is k- 
times that of a single random walk [Eq. (0}]: A qs (Sfc) = 
k. Using Eqs. © and {BJ we find: 

4fc 2 - 1 

A{S k ) = (Al) 

which correctly reproduces A(5 a ) of Eq. JT^} for k = 2. 
In the case of four mutually avoiding walks 

A(S 4 ) = f (A2) 

Notice that the vertex S c of Fig. 0] is also formed by four 
random walks, however not all mutually avoiding and its 
bulk conformal dimension is A(S C ) — 35/24 < A(S , 4 ). 
A(S C ) is smaller as more configurations are available for 
the star copolymer when partial overlapping between 
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walks is allowed, as in S c . Had we put no restrictions ated. It should be emphasized that for two dimensional 
on the order of the loops in the construction of the star copolymers the order of the constituing walks does 
model both vertices S c and S4 would have been gener- matter. 
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